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Abstract. The Wroblewski parameter is a convenient indicator of strangeness production and can be 
employed to monitor a signal of quark-gluon plasma production : enhancement of strangeness production. 
It has been shown to be about a factor two higher in heavy ion collisions than in hadronic collisions. Using 
a method proposed by us earlier, we obtained lattice QCD results for the Wroblewski parameter from our 
simulations of QCD with two light quarks both below and above the chiral transition. Our first principles 
based and parameter free result compare well with the A- A data from SPS and RHIC. 

PACS. 12. 38. Mil Quark-gluon plasma - 12.38.Gc Lattice QCD calculations 



1 Introduction 

As with many signals of quark-gluon plasma (QGP) pro- 
duction in relativistic heavy ion collisions, the basic idea 
behind enhancement of strangeness production jlj as a 
QGP signal is very simple. Recognising the fact that the 
strange quark mass is smaller than the expected transi- 
tion temperature whereas the mass of the lowest strange 
hadron is significantly larger, it was argued that the pro- 
duction rate for strangeness in the QGP phase, o-qgp(ss) 
is greater than that in the hadron gas phase, <jhg{ss). 
While this energy threshold argument for strangeness pro- 
duction in the two phases is qualitatively appealing, one 
has to face quantitative questions of details for any mean- 
ingful comparison with the data. Applications of pertur- 
bative QCD needs a large scale which could be either the 
temperature of QGP or the mass of the produced strange 
quark-antiquark pair. Since the temperature of the plasma 
produced in RHIC, or even LHC, may not be sufficiently 
high for perturbative QCD to be applicable and since 
the strange quark mass is also rather low , estimation 
of strangeness production by lowest order processes like 
gg — > ss could be misleading. Indeed, it is now well-known 
that even for the charm production, the next order cor- 
rection to gg — > cc is as large as the leading order; such an 
order by order computational approach may be hopelessly 
futile for the lot lighter strange quark. 

A variety of aspects of the strangeness enhancement 
have been studied and many different variations have been 
proposed. One very useful way of looking for strangeness 
enhancement is the Wroblewski parameter Defined as 



the ratio of newly created strange quarks to light quarks, 



2(ss) 



dd) 



(1) 



the Wroblewski parameter has been estimated for many 
processes using a hadron gas fireball model 3 . An in- 
teresting finding from these analyses is that A s is around 
0.2 in most processes, including proton-proton scattering, 
but is about a factor of two higher in heavy ion collisions. 
An obvious question one can ask is whether this rise by 
a factor of two can be attributed to the strangeness en- 
hancement due to quark gluon plasma and if yes, whether 
this can be quantitatively demonstrated by explicitly eval- 
uating the Wroblewski parameter in both phases. Alter- 
natively, one could just study how different the prediction 
actually is and learn about other physics effects from its 
comparison with data. We show below how quark number 
susceptibilities, obtained from simulations of lattice QCD, 
may be useful in answering such questions. Since these 
simulations correspond to equilibrium situations, one needs 
certain extra assumptions which we also discuss briefly 



2 A s from Quark Number Susceptibilities 

Quark number susceptibilities (QNS) can be calculated 
from first principles using the lattice formulation. Assum- 
ing three flavours, u, d, and s quarks, and denoting by [if 
the corresponding chemical potentials, the QCD partition 
function is 
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Fig. 2. Typical continuum extrapolation results. 



Note that the quark mass and the corresponding chemi- 
cal potential enter only through the Dirac matrix M for 
each flavour. We use staggered fermions and the usual 
fourth root trick |3] to define M for each flavour. Defin- 
ing no = Hu + £id + Ms and ^ 3 = /Lt u - fid, the baryon and 
isospin densities and the corresponding susceptibilities can 
be obtained as: 



Td\nZ 
V dfM ' 



Xij 



T d 2 lnZ 
V dfiidfij 



(3) 



on the methods and/or models which aim to explain the 
large deviations of the lattice results for pressure P(/x=0) 
from the corresponding perturbative expansion. Here we 
will be concerned with the Wroblewski parameter which 
we |H] have argued can be estimated from the quark num- 
ber susceptibilities: 



QNS in eq. JJJ) are crucial for many quark-gluon plasma 
signatures which are based on fluctuations in globally con- 
served quantities such as baryon number or electric charge. 
Theoretically, they serve as an important independent check 



X, 



2Xs 



Xv 



Xd 



(4) 



Note that the lattice simulations yield real quark num- 
ber susceptibility whereas for particle production its imag- 
inary counterpart is needed. Indeed, \ s above too needs 
the latter. However, one can relate the two and thus justify 
the use of lattice results in obtaining \ s . Briefly, the argu- 
ment jH| is as follows. Fluctuations in physical quantities, 
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Fig. 3. Quenched QCD results in continuum limit. 
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Fig. 4. Comparison of the corresponding A s with RHIC and SPS experiments. 
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Fig. 5. A s as a function of temperature for full and quenched QCD. 
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described by a perturbation in time, can be related to a 
generalized susceptibility for the corresponding operator 
for it. This is complex in general. Its imaginary part can 
be shown to determine the dissipation, i.e., the production 
of strange quark-antiquark pair in our case. From the gen- 
eral properties of these susceptibilities, a Kramers-Kronig 
type relation between their real and imaginary parts can 
be obtained. Finally, making a relaxation time approxima- 
tion {lot 3> 1), one finds that the ratio of the imaginary 
parts is the same as that of the real parts. 

In order to use eq. to obtain an estimate for com- 
parison with experiments, one needs to compute the corre- 
sponding quark number susceptibilities on the lattice first 
and then take the continuum limit. All susceptibilities can 
be written as traces of products of the quark propagator, 
M (rrig), and various derivatives of M with respect to 
\x. With m u — md, diagonal Xu' s can be written as 



Xo = ^[{Oi{m u ) + \o n {m u ))} 
X3 = ^ (0 2 {m u )) 



(5) 
(6) 
(7) 



Here 2 = Tr M^M^-Tr M~ 1 M' u M ~ 1 M' u , and O n {m u ] 
(Tr M~ 1 M' V )' 2 . The traces are estimated by a stochas- 



tic method: Tr A = Y^ZiR\ARi/2N v , and (Tr A) 2 = 

2 Ef>j=i( Tr ^M Tr -4) J /i(i - 1), where R, is a complex 
vector from a set of N v , subdivided further in L indepen- 
dent sets. We use typically N v = 50-100. 

Figure 1 displays results [Hj for the susceptibilities as 
a function of temperature in units of T c , where T c is the 



transition temperature. Normalized to the corresponding 
ideal gas results on the same lattice, i.e, the infinite tem- 
perature limit, results for QCD with two light dynamical 
quarks of mass 0.1 T c are shown as points whereas the con- 
tinuous curves correspond to the results in the quenched 
approximation. Note that the latter amounts to dropping 
the fermion determinant term in simulations which be- 
come orders of magnitude faster, and hence more precise, 
than the full QCD simulations. The valence quark mass 
m v , appearing in eqs. ©-(01 is shown in the figure in 
units of T c . Note that T c in these two cases differ by a 
factor of 1.6 but the results for the corresponding dimen- 
sionless susceptibilities as a function of the dimensionless 
ratio T/T c differ by a few per cent only. Such a mild 
dependence on the number of dynamical flavours in the 
thermodynamic quantities has been a known feature in 
the temperature region away from the transition. Indeed, 
since the nature of the transition does depend strongly on 
the number of dynamical flavours, one expects significant 
differences near T c . Encouraged by this behaviour, we in- 
vestigated the continuum limit for the quenched case by 
increasing the temporal lattice size from 4 to 14 in steps 
of two and extrapolating to infinite temporal lattices. The 
spatial lattices were also increased to maintain the as- 
pect ratio constant. Figure 2 shows typical results of such 
continuum extrapolation at T = 2 T c . The continuum re- 
sults for the light quark susceptibility thus obtained in 
the quenched approximation are exhibited in Figure 3 for 
small m v . The bands marked by 'HTL' and 'NL' show 
the analytic results of |7j obtained in successively better 
approximations respectively. 

The strange quark susceptibility in the continuum limit 
was obtained from the same simulations by simply choos- 
ing m v /T c — 1 (in both full and quenched QCD in view 
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of Figure 1. Using eq. (@J), A S (T) can then be easily ob- 
tained. These were extrapolated to T c by employing simple 
ansatze. The resultant A S (T C ) in quenched QCD is shown 
in Figure 4 along with the results obtained from the anal- 
ysis of the RHIC and SPS data in the fireball model[3]. 
The systematic error coming from extrapolation is shown 
by the brackets. The agreement of the lattice results with 
those from RHIC and SPS is indeed very impressive. 

The nice agreement needs to be treated cautiously, 
however, in view of the various approximations made. Let 
us list them in order of severity. 

— The result is based on quenched QCD simulations and 
extrapolation to T c . As seen from Figure 1, the quark 
number susceptibilities, and hence A S (T), are expected 
to change by only a few per cent. Since the nature of 
the phase transition does depend strongly on the num- 
ber of dynamical quarks, a direct computation near T c 
for full QCD is desirable. We are currently making such 
a computation and have some preliminary results for 
full QCD with two light dynamical quarks for lattices 
with four sites in temporal direction. These are shown 
in Figure 5 along with the continuum quenched results 
for X S (T) and the band for experimental results. While 
the emerging trend is encouraging, further exploration 
with varying strange quark mass, temporal lattice size 
(to obtain continuum results) and spatial volume is 
still necessary. 

— The experiments at RHIC and SPS have nonzero albeit 
small fj, whereas the above result used /i — 0. Based 
on both lattice QCD and fireball model considerations, 
A s is expected to change very slowly for small \i. This 
can, and should, be checked by direct simulations. 

— As argued above, particle production needs the imag- 
inary counterpart of what one obtains from simula- 
tions. The relation between the ratios of real and imag- 
inary parts was obtained under the assumption that 
the characteristic time scale of quark-gluon plasma are 
far from the energy scales of strange or light quark pro- 
duction. Observation of spikes in photon production 
may falsify this assumption. 
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3 Summary 



Quark number susceptibilities contribute in many differ- 
ent ways to the physics of the signals of quark-gluon plasma 
in heavy ion collisions at SPS and RHIC. They can be 
obtained from first principles using lattice QCD. This of- 
fers a quantitative control and check of these signals and 
thus QGP itself. In particular, the continuum limit of Xu 
and Xs, which we obtained in quenched QCD, leads to a 
temperature dependent Wroblewski parameter, A S (T). Its 
extrapolation to T c appears to be in good agreement with 
results from SPS and RHIC. First full QCD results near 
T c confirm this as well, although many technical issues, 
e.g, finite lattice cut-off or strange quark mass, need to be 
sorted out still. 



